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ABSTRACT 


The  modified  hypernetted  chain  theory  (MHNC)  for  density  profiles 
of  fluids  in  contact  with  a  hard,  smooth  surface  is  found  to  posess  a 
remarkably  simple  and  interesting  property:  The  accuracy  of  the  MHNC  and 
the  universality  of  the  bridge  functions  for  the  density  profiles  can 
be  tested  without  resort  to  any  detailed  solution  of  the  integral 
equations.lt  is  given  by  the  degree  of  universality  of  the  bulk  bridge 
parameter  (e.g.  rj  for  the  Percus -Yevick  hard-sphere  bridge  function) 
when  expressed  in  terms  of  the  second  partial  derivative  of  the  free 
energy  with  respect  to  the  density 

F  --pZ/2  ^  f(j),T)/  dp  | 

The  function  rj  (F)  as  obtained  from  MHNC  calculations  for  bulk 
simple  fluids  is  found  to  be  remarkably  independent  of  the  pair 
potential. As  an  example  for  the  accuracy  of  the  resulting  method,  the 
one  component  plasma  near  a  hard  wall  is  discussed. 


The  pair  correlation  functions  and  the  equation  of  state  of  bulk 
are  given  with  a  high  degree  of  accuracy  by  the  MHNC, which  is  based  on 
the  universality  of  the  short  ranged  part  of  the  bridge  functions  [1], 
However,  a  similarly  accurate  theory  for  the  inhmogeneous  fluid  is  not 
yet  available. In  particular  this  is  the  case  of  the  interface  of  a  ionic 
solution  near  a  charged, smooth  hard  wall  [2], This  is  a  simple , highly 
idealized  model  of  the  interface  of  a  ionic  fluid  and  a  charged  solid 
but  it  has  been  extensively  used  in  electrochemistry  [ 3 ] .Recently , as 
part  of  an  extensive  study  of  the  one  component  plasma  (OCP)  [4], the 
interfacial  properties  of  the  OCP  and  hard  surfaces  [$-8]  were  studied. 

The  comparison  between  the  accuracy  of  the  same  approximation  in 
the  bulk  or  at  the  interface  shows  that  in  general, a  good  approximation 
for  the  bulk  phase  is  not  as  succesful  for  the  interface . For  example, 
the  HNC  is  an  excellent  theory  for  the  bulk  electrolytes  and  OCP  [1,4] 
However  it  is  a  much  poorer  theory  for  the  charged  interface .The  reason 
is  that  it  underestimates  the  density  oscillations  caused  by  the 
excluded  volume  effects. The  error  produced  by  this  underestimation  is 
much  more  serious  for  the  flat  interface  case. The  inclusion  of  bridge 
diagrams  corrects  in  part  this  defficiency,  and  should  yield  a  much  more 
accurate  theory  for  the  interface. 

As  a  first  step  in  the  systematic  improvement  of  the  flat  wall  HNC 
similar  to  what  has  been  done  in  the  bulk  case  [1] ,  the 
present  paper  provides  a  proof  of  the  universality  of  the  short  ranged 
part  of  the  bridge  functions  for  the  density  profiles  near  a  flat, 
smooth  interface .The  analysis  of  the  MHNC  scheme, and  the  availability  of 
MHNC  results  for  the  bulk  OCP  enables  us  to  establish  this  result 
without  having  to  particularize  to  any  specific  density  profile  in 


detail. For  the  restricted  (dielectric  continuum)  primitive  model  of 


£ 


K 


ionic  solutions  a  similar  analysis  is  also  possible. 

A  fluid  in  contact  with  an  impenetrable  wall  can  be  considered  [9] 
as  a  limiting  case  of  a  mixture  in  which  one  of  the  components  (which  we 
label  w)  grows  in  size  but  dwindles  in  concentration,  so  that  the 
properties  of  the  bulk  phase  of  the  fluid  remain  unchanged. Let  R  be  the 
hard  core  radius  of  the  w  particle  (the  "wall"),  with  density  p  .Then 

*  W 


the  limit  is  [9,10' 


-»0 


R  ~>qo  and  «  R  -5*0 . 

w  )  w  v/ 


We  use  the  Ornstein  Zernike  (OZ)  equation  as  our  starting  point. We 
show  that  in  the  planar  limit  the  bridge  functions  of  the  MHNC 
approximation  can  be  determined  from  a  universal  function. The  criterion 
to  determine  the  parameters  of  this  function  are  given  in  section  II. 

The  remarkable  (but  perhaps  not  unexpected)  fact  is  that  nearly  all  the 
bridge  function  parameters  are  already  determined  by  the  solution  of  the 
bulk  MHNC  equations .The  only  free  parameter  turns  out  to  be  a  shift 
parameter,  ^  .which  determines  the  position  of  the  wall  in  the  reference 
system  relative  to  that  of  the  system  under  consideration.  This  special 
property  of  the  bridge  functions .when  analyzed  in  the  context  of  the  sum 
rules  for  the  density  profiles  (discussed  in  section  III),  leads  to  a 
remarkably  simple  and  interesting  result  :The  accuracy  of  the  MHNC 
theory  and  the  universality  of  the  bridge  functions  for  the  density 
profiles, can  be  tested  without  having  to  solve  the  equations.lt  is  given 
by  the  degree  of  universality  (  i.e.  the  independence  of  the  pair 
potential)  of  the  bulk  bridge  parameters  (for  example,  ij  of  the 
Percus-Yevick  hard  sphere  functions),  when  expressed  in  terms  of  the 
second  partial  derivative  of  the  excess  free  energy  with  respect  to  r j. 


•VJW.VJV w. 
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This  test  is  psi'UL'i.cu  in  section  i.  .  , inere  tne  univeisaiitv  is  snown  zo 

be  very  accurate . (see  figure  1) .As  an  example  of  the  resulting  method 
and  its  accuracy  we  consider  in  section  V  the  density  profile  of  the  OCP 
near  a  flat  wall. The  Percus  Yevick  hard  sphere  bridge  functions  for  the 
density  profiles  are  briefly  discussed  in  section  VI  while  some 
implications  of  the  present  work  are  mentioned  in  the  concluding  section 


VII. 


II -Ornstein-Zernike  Equation  and  MHNC  Closures  in  the  Planar  Limit 


Consider  a  mixture  of  particles  labelled  l,2,..m,w  .with  total 
density 


*  I 


r2."‘ 


.  /V-N/V, 


L-  t 


-l 


temperature  Z-(kftT)  ,and  mole  fractions  x.-N.  /N  ,x  -N ^/N.In  the  limit 
x  ->0  the  OZ  equation  decouples  into  two  parts  [9]: the  bulk  part 

W 


h.,(r)-c.  (r)  +p  7  x  (  dr'h  (|r-r*  |)c  (r') 

fc>  ~  ~  L)  ~  i*  -  -  el¬ 

and  the  wall  part  *•"' 


(2.1) 


h  (r)-c :  (r)  +p  )  x  f  dr'h  ( |  r-r*  |  )c.  (r ' ) 
w;  ■»  w;  ~  )  /  f  J  —  wji  -  —  ~ 


(2.2) 


where  h.(r)  is  the  total  correlation  function,  and  c  (r)  is  the  direct 
'i  -  ~ 

correlation  function  for  the  bulk  fluid,  h  (r)  and  c  (r)  are  the 
correlations  involving  the  w  particle . Consider  now  the  planar  limit:  let 
r-R  +x  (eventually  R  -»o«)  and  denote 

vf  w 


h  ,(r)  yh.(r) 


c  (r)  -»c  (r) 


Then  [9] 


oO 


h^.(x)«c  -(x)  +p  ^Xj2ft  Jdt  h^(t)jCs  sc^,(s) 

I  i.M 


(2.2') 

JU 

u-tt 

where  we  are  implicitly  using  hard  core  exclusion  conditions  for  the 
wall 


g.(x)-h;(x)  +1-0, 


x<R 


(2.3) 


where  R-is  the  radius  of  the  hard  core  of  i. 

I 

The  fluid  now  occupies  Che  right  half  space , x>0 , and  jox^ g  •  (x)  gives 
the  probability  of  finding  the  center  of  particle  i  at  a  distance  x  from 
the  surface  of  the  wall.  The  function  g.(x)  represents,  therefore,  the 
density  profile. When  the  wall  is  charged  and  the  system  has  ions,  then 
(2.2')  has  to  be  modified  to  avoid  divergencies  produced  by  the  Coulomb 
potential  [11] 


oO 


Jds 


sc  (s) 

U 


oO 


Jds 


sc  (s)  +  SQ .  Q( | x- t I +x+t) 
4L  1  * 


(2.2") 


where 


St 

sc^(s)  -  sc^  (s)+  SQ^Qj^/s  (2.4  ) 

The  charged  system  must  satisfy  the  bulk  electroneutrality 
condition 

Z*;Qc°  (2-5> 

b  “  1 

and  the  wall  electroneutrality  condition 


jo^x  .Q-Jdx  g.(x)  -  -E0/4TT 


(2.6) 


where  E  is  the  electric  field  at  the  surface  of  the  wall, 
o 

The  exact  diagrammatic  expansion  of  the  fluid  pair  correlation 
functions  supplements  the  OZ  equation  with  the  closures  [1,4] 


The  hvpernetted  chain  approximation  HNC  is  obtained  by  setting  B..i,r;=U. 

Note  also  that  (2.7)  is  an  HNC  for  the  effective  potential 
A  (r)+B  (r)/£ . Presently ,  the  MHNC  is  the  most  accurate  theory  available 
for  the  pair  correlation  functions  of  the  f luld[ 1 , 12 - 16 ] . The  MHNC  scheme 
is  based  on  the  short  range  universality  of  the  bridge  functions: In 
practice,  one  considers  a  reference  system  (Superscript  'o')  of  hard 
spheres  of  radii  R°,R°,  with  the  same  relative  concentrations  x  x  , 

i  w  w- 

temperature  &  and  density  jo  as  the  original  system. The  basic  MHNC 

approximation  is 


B  (r)  »B  (r ,R° _ R° ,R°) 


(2.8) 


where  the  reference  hard  core  radii  R?,R°  are  adjustable  parameters  to 

l  W 

be  determined  by  imposing  consistency  criteria. In  particular,  the 
following  MHNC  ' virial-energy'  consistency  requirement  can  be  imposed 
locally  [13-16]  (without  having  to  consider  different  temperatures  or 
densities)as  a  set  of  m+1  equations  for  the  m+1  unknown  radii  R*(S,o): 


/  x.  x.  fdr  [  g.(r)  -g°(r)]  <*B;,(r)/6R^  -0 

^ - t  ‘  ij  ~  *j-  ‘p  H 

j  ^  :  *• "  wv'i'*r  k-l,..ra,w 


(2.9) 


The  HNC, with  B..  -0,  is  the  member  of  the  lowest  order  of  this  class  of 
approximations  satisfying  ’virial-energy'consistency  [13,14,17], 

Detailed  results  [1,12-16,18-21]  show  that  the  MHNC  reproduces  the 
computer  simulation  results  for  a  large  variety  of  fluids  within  the 
statistical  error  of  the  simulations .The  results  for  multicomponent 
classical  plasmas  are  nearly  as  accurate  [12], while  somewhat  less 
accurate  results  are  obtained  for  electroytes  and  molten  salts  [20]. 

In  any  event,  the  MHNC  results  in  a  significant  improvement  of  the  HNC 
theory. The  accuracy  of  the  local  MHNC  theory  based  on  the  universality 


perturbation  theory  for  the  equation  of  state,  were  recently  studied  in 
detail  [13-14] . 

The  interesting  feature  of  the  planar  limit  is  that  the  closures 
(2.7)  decouple  from  the  bulk  part  in  the  same  way  as  the  OZ  equation: 


log  (1+h  .  (r) )-  -  [iU..(r)+B.  (r)  ]  +  [h.  (r)  -c..(r)  ] 

"J  "1  >•) 

i ,  j-1 . .  . m  (2.10) 

and  for  the  wall  part 


log  (1+h  •  (x) )-- (x)+B • (x) ]+[h . (x) -c- (x) ] 

v  /  v.  L  tv. 

i-l...m  (2.11) 

where  ^.( r )  -><}.  (x)  ,B^(r)  -*>B;  (x)  . 

When  employing  the  MHNC  scheme  for  the  bulk  part.eqs.  (2.8,9) 
remain  unchanged,  except  that  the  w  component  is  not  taken  into 
account. In  other  words, we  get  the  usual  MHNC  bridge  functions  for  the 
bulk  phase  with  the  m  bridge  parameters  R°  (11,^)  .Taking  the  planar  limit 
in  the  reference  system  (using  now  r-R^.  +z-Rw+)v+x)  ,  we  obtain  the 
reference  bridge  functions 

g*(r)-*g‘.(z)  >g’(x+>,R’  (S,p) . R°  (JS.p)) 

B^r)->BQ(z)  — B?  (x+^, ,  R°  (fi.p) . R^(fi,j>))  (2.12) 

so  that  the  (m+l)th  equation  in  (2.9)  becomes 


from  which  we  find  the  remaining  bridge  parameter  (the  shift  parameter) 


III-MHNC  Equations  tor  cne  Density  Protnes  and  Bridge  rune  cions :  Sum 

Rules 


Combining  (2.2")  and  (2.11)  we  arrive  at  the  following  equations 


for  the  density  profiles: 


WN 

y  /  /• 

log[l+  h  -  (x) ]-  p  J  jdt  h^(t)jds  sck(s)  -ilQJjT  (x)  -  &  <|  .(x) 


-  b^(x)  , x^R^ 
(3.1) 


[1+  h.(x)]-0,  x<Rt 


where  lj(x)  is  the  electrostatic  potential  in  the  half  space  X^C 


\jf(x)-4rff>  Jdt  (x-t)  ^Jx^h^ (t) 


which  satisfies  the  boundary  conditions 


<ty/^x  j  --Ee  ;  l|J(x-^oo)-0 


(3.2) 


(3.3) 


i.(x)  is  the  non-electrostatic  , non-hard  core  part  of  the  interaction 
1  v 

between  particle  i  and  the  wall.Eqs  (3.1-3)  also  apply  to  the  situation 
when  the  system  is  immersed  in  a  uniform  neutralizing  background  of 
charge  density  p^  :  It  this  case,  however,  (2.5  )  and  (2.6)  have  to  be 


modified 


f  Lx  Q  "K 

h  S  I 

f>  7*.  Q  Jdx  h-(x)  -  -Eo/41T 


(3.4) 


(3.5) 


When  che  bridge  functions  of  a  reference  system  are  employed,  then  using 
the  property 


bB°t(x+X  )/*>  -bB’(x+X  )/6x 


we  arrive  at  the  following  condition  for  the  shift  parameter^: 


W' 


^  xvj(ix  [  g.(x)  -g°(x+),  )]  'bB?(x+^  )/bx  >0 
c-'  i-1,  .  .m  (3.6) 


It  should  be  noted  that  the  HNC  approximation  for  the  density 
profiles  (i.e.  B^(x)—0  without  any  additional  statement  about  the  bulk 
B.^(r)'s)  is  obtained  also  from  Wertheim's  equation  [22]  :  For  a  one 


component  system 


VV08f(£}  +k£  > -/>]  (3-7> 


under  the  assumption 


c(r , r ' )-  c  ( | r-r ' | ) 


(3.8) 


Employing  the  method  of  Carnie  et  al  [23],  we  derive  from  (3.1)  and 
(3.5)  the  following  sum  rule 


V* 


^  xcg,(V  &Ec/8rr^  +i/2[i+&(&p/dp)T] 


Lk 


(.  • ' 


*<*  oO 


+  ^xiJdx  g^(x)  fi^.(x)/6x  +^x-jdx  g^(x)  <)B^(x)/bx 


(3.9) 


V  -  I 


If  the  exact  bulk  direct  correlation  functions  (dcf)  and  the  exact 


bridge  functions  B;(x)  are  employed,  then  eq.(3.9)  is  exact.lt  should 
be  compared  to  the  following  exact  relation  that  does  not  involve  the 
bridge  functions  explicitly  [2,24,25] 
j  xiS.(R.)  +&  (Jx^)  <£(x-0). 


rr  5*2. 

+  J  *vJdx  g^(x)  11  "bji -(x)/dx 


t  *  i 


(3.10) 


comparing  these  equations  we  find  the  following  bridge  function  sum  rule 
for  the  density  profile 


v<« 

x-Jdx  g.(x)  l>B-(x)/dx-  ^-l/2[l+B(^P/bp)T]  uL< 

(3.11) 

When  approximate  bridge  functions  e.g  B°t(x)  and  /or  approximate 
bulk  dcf  (c  (r))  are  used,  ther^eq.  (3.9)  is  still  valid  provided  that 
the  bulk  inverse  compressibility 

r 


is  taken  from  the  compressibility  equation  of  state  which  corresponds  to 
the  bulk  dcf's  that  are  employed. When  B^(x)-0  then  we  recover  the  more 
familiar  HNC  sura  rule  [23], In  the  case  of  the  mean  spherical 
approximation  (MSA)  the  closure  for  the  dcf  is 


c  ,(x)  -•  A.(x)  ,  x^R; 

•  w 


(3.12) 


for  which  we  obtain  from  (2.11) 


JM5 


(3.13) 
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g  (x)  'bB  (x)/c)x|  -g  .  (x)dg-(x)/dx  -bg(x)/bx 
v  1  1 


and  from  (3.9) 

n*'  »V>  J 

1/2  +<i  (ZX‘Qc)  ftEe/8lTf  +1/2[R(Sp/3o)t 
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x(Jdx  S ( x )  &  3<|.(x)/3x 


(3.14) 


Despite  the  extensive  work  on  the  MSA  for  density  profiles 
[2,11,26],  this  general  MSA. sum  rule  has  not,  to  our  knowledge  been 
published  in  the  literature . This  MSA  sum  rule,  as  well  as  the  general 
sum  rule  (3.9),  are  very  useful  in  understanding  the  results  of  the 
different  approximations. 

Returning  to  the  MHNC  approximation  with  the  hard  sphere  reference 
bridge  functions  B! (x+X) ,  we  write,  from  (3.6) 


xv-jdx  g.(x)  bB?(x+)0/»x  -  ^xjjdx  g»(x)  bB ?(x)/bx 

wt  3 

x.jdx  g!(x)  ^B^xJ/bx 


i=.  o 


(3.15) 


Using  (3.11)  applied  to  the  reference  system,  and  also  (3.9),  we 
obtain  the  following  MHNC  result  for  the  density  profiles: 


^  xtgt(Rt)  (2,x-Qc)  £(x-0)-  BE,/8fT/> 

- 1  ;i1 

Ox.  dx  g.(x)  ft  'hi.  (x)/Sx  +l/2[l+fi(&P/dp)  ]fc>jlK  +[ft(P/^) 

T-1  J  t  ’  *■  ^  NkMOTC 

vW  ^ 

-l/2[l+B(3P/dp^  ^  -Jx./dx  g!  (x)  ft  ^(x)/dx 


(3.16) 


where  ,  for  all  practical  purposes  we  may  ignore  the  last  term  (i.e.^R^ 
of  the  right  hand  side  [27], 

Comparing  (3.16)  with  (3.9)  and  (3.10),  we  find  that  in  all 
cases  when  the  bulk  dcf's  and  equation  of  state  are  given  correctly  by 
the  MHNC  (as  it  is  the  case  for  simple,  one  component  fluids), a  criteria 
for  the  performance  of  the  accuracy  of  the  MHNC  for  the  density  profile 
of  the  inhomogeneous  case  is  provided  by  the  validity  of  the  following 
bulk  relation  (the  'bridge  sum  rule') 


[B(p/j))]-i/2[i+&(ap/ay)T]  -  [&(p/j>)r-  i/2[i+s(3p/dj>)T 


(3.17) 


Denoting  by  f(p,T)-ftF  /N  the  excess  free  energy  per  particle 
then  we  may  re-write  (3.17)  as 


-j)V2  9  f(y,T)/dp4  |—^/2  'dV(j>,T)/3p4!_r 

(3.18) 


Thus,  an  optimized  MHNC  map  ( i  .  e  .  R”  (6 , jj) )  of  the  system  on  the 
reference  system,  such  that  the  optimized  bridge  parameters  obtained 


from  (2.9)  reproduce  accurately  the  structure  functions  of  the  system, 

will  perform  as  well  for  the  density  profiles,  provided  that  the 
optimized  reference  bulk  system  and  the  given  bulk  system  have  the  same 
second  partial  density  derivatives . This  result  is  an  unexpected  type  of 
correlation  between  the  bulk  bridge  parameters. 


IV-Test  of  the  Bridge  Sum  Rule : Universal i ty  of  the  Wall  -  Particle  3ridge 

Functions 


For  bulk  fluids  in  three  dimensions  the  exact  bridge  functions  are 

not  known: However  the  universality  hypothesis  can  be  tested  very 

effectively : The  bridge  functions  obtained  from  the  Percus-Yevick  theory 

for  hard  spheres  reproduces  the  structural  correlation  functions 

obtained  by  computer  simulation  of  a  large  variety  of  cases,  with  quite 

disparate  interaction  potentials,  from  hard  spheres  to  charged  plasmas 

[ 1 , 12- 21 ] . The  corresponding  MHNC  calculations  employ  the  PY-hard  sphere 

bridge  functions  with  a  single  parameter  rj  (which  is  the  hard  sphere 

packing  fraction  in  the  PY  approximation) .As  a  result,  we  have  available 

the  values  of  rj  (^>,T)  for  the  set  of  potentials  for  which  detailed  MHNC 

calculations  have  been  performed. These  include  the  hard  spheres, 

-  w  -  <•  _  \ 

Lennard-Jones  (LJ  12-6)  ,r  ,r  and  r  (the  OCP)  potentials [ 1 , 6 , 18 , 19 ] . 
As  already  discussed  in  detail  elsewhere  [13]  when  the  bridge  parameter 
is  expresed  in  terms  of  the  entropy  S~-S  /Nk  , then  ^(S)  forms  a 
universal  line,  independent  of  the  potential  (fig  1  of  ref".  13). This 
represents  a  graphic  demonstration  of  the  accuracy  of  the  variational 
thermodynamic  perturbation  theory. The  result  (3.17)  has  not  been 
anticipated. 

The  test  of  the  validity  of  (3.17)  and  of  the  universality  of  the 
bridge  functions  for  the  density  profiles  can  be  performed  together  by 
examining  the  universality  of  the  bulk  bridge  parameter  n  ,when 
expressed  in  terms  of 

F—  [  & ( P/j> )  ]  - 1/2  [  1+&  ( S P/dp  )t  ] 


v _ 


& 


/2  3  f(p,T)/dp 


i-».  I) 


for  different  potentials.  It  is  to  be  emphasized  that  this  broad 


interpretation  of  the  universality  r\_  (F)  requires  only  that  a  set  of 


bridge  functions  from  which  this  relation  is  obtained  does  accurately 


reproduce  the  bulk  pair  correlation  functions. In  figure  1  we  plot  the 


results  for  i\  (F)  for  different  potentials  in  the  entire  fluid  range, 


r>£<0 .45 .  Indeed,  the  universality  is  satisfied  to  a  ramarkable  degree. 


Better,  in  fact  than  that  of  ry[S),The  hard  sphere  line  in  figure  1  is 


actually  obtained  from  simple  analytic  expressions.  For  hard 


spheres , using  the  packing  fraction  ^-T7/6  pd*  (where  d  is  the  hard 


sphere  diameter),  we  use  the  Carnahan-Starling  equation  of  state  [29]  to 


(4.2) 


This  expression  provides  also  an  excellent  fit  for  the  MHNC 


calculations  of  Tsai  [18], who  obtained  the  representatior^for  the 


parameter 


ti-  (0.982913-0.02271^  +0.02449  }  ) 


(4.3) 


The  universal  behaviour  of  iq(S)  and  n^F)  for  simple  one  component 


fluids  suggests  a  new  interesting  possibility  to  correlate  the  bulk 


properties  of  simple  classical  fluids  using  a  relation  of  the  type 


S(F)-  universal  function 


(4.4) 


V-An  Example: The  one  Component  Plasma  in  Contact  with  an  Impenetrable 

Wall 


Monte  Carlo  (MC)  simulations  of  the  surface  properties  of  the  OCP 
were  performed  by  Badiali  et  al  [5], These  authors  present  the  results 
for  an  near  an  impenetrable  hard  wall  for  five  values  of  the  plasma 
coupling  parameter 

2.  .  'll 

P-J* Q  /a*  -1,10,20,30.70  a  -[3/4TTf  ] 

In  addition  to  the  MC  result  they  also  present  results  for  the  HNC 
closure  for  the  density  profile  using  the  MC  results  for  the  bulk  dcf, 
and  results  for  the  MSA  for  the  density  profile  [26]  using  the  soft  MSA 
[30,17]  bulk  dcf's.We  will  refer  to  those  results  as  MC/MC.HNC/MC  and 
MSA/SMSA . 

s* 

For  the  case  at  hand,  (x)-0 , Eq-0 , R-0 ,  and  in  the  notations  of  ref 
[5],  &QU(0)--fize*^  .Denoting 

-a  Y”  -g(0)  in  general  and  4  Eze*^  -gZ(0)/2  for  the 

MSA, we  obtain  the  following  results  for  the  different  theories 


o  y  -l/2[l+E(dP/dp)T]0v)t< 

r*c 

*  ^  -l/2[l+S<6P/dp)T]^ 

Sms* 

-l/2[l+fi(dP/0p)  ]^tH 

Mn-c. 

M  -fi(P/p)  ] 

*2  -&(p/p)  -(F-F0) 


for  HNC/MC 


for  MSA/SMSA 


for  HNC/HNC 


for  MC/MC 


for  MHNC/MHNC 


(5.1) 


Recall,  that,  essentially,  MHNC/MHNC^MHNC/MC  and  that  the  OCP  results 
give  ^  ( r1 -l)-0 . 1 ,  ( r*  ~70)-0 . 38  (see  fig.l  in  ref  6},Using  fig.l  with 

^F«F-F°  given  in  relation  to  the  hard  sphere  line  we  find 

for  V  -  1  ,  10  ,  20  ,  30  ,70 

we  find^F  -  0  ,  0.20,  0.25  ,  0.15  ,  -0.25  respectively .  Recall  also 

that  for  the  bulk,HNC~MC  for  r-1,  while  [17,31],  HNC-SMSA  for  f>10. 

The  reults  of  the  different  theories  for  ^  ^1  are  given  in  table  l.Note 

the  high  accuracy  obtained  by  the  MHNC/MHNC  when  compared  to  the  MC/MC. 

Note  also  that  without  succeeding  to  reproduce  the  pronounced  structures 
of  the  MC/MC  density  profiles, all  the  theories  of  the  type  X/HNC  or 
X/SMSA  with  X-HNC  or  MSA  will  give,  for  large  plasma  parameter  TP  ,  much 
better  results  for  the  total  potential  drop  &  ^  than  the  corresponding 
X/MC.This  is  a  negative  feedback  artifact  of  the  thermodynamic 
inconsistency  of  the  HNC  and  SMSA  results  for  the  bulk  OCP,  by  which 
they  give 

l/2£  (dP/dp 0.3  Hwich  happens  to  agree  with  &  (P/p)  -  0.3  I-* 

This  situation  is  very  similar  to  what  happens  in  the  bulk  OCP  where 
both  the  HNC  and  the  SMSA  give  very  good  results  for  the  potential 
energy  despite  )  the  fact  that  the  corresponding  pair  correlation 

functions  are  far  from  reproducing  the  MC  results.  The  MHNC/MHNC  results 
satisfy  a  stringent  moment  test,  i.e.(3.6)  which  as  in  the  bulk  case 
also  assures  the  validityof  the  *  y  (<->  potential  energy)  test  in  table 


VI -Construction  of  the  Bridge  Functions  for  the  Density  Profiles 


Despite  the  fact  that  we  have  used  the  Pi  hard  sphere 
representation  for  the  bulk  bridge  functions  in  constructing  the  MHNC 
scheme  for  the  inhomogeneous  problem, the  universality  criteria  for  the 
bridge  functions  for  the  density  profiles  in  the  inhomogeneous  fluid  and 
eq.  (3.17)  are  valid  in  all  cases.  The  universality  of  n(F)  tells  us 
that  there  is  a  similarly  universal  representation  of  tne  bridge 
functions  for  the  density  profiles,  but  it  does  not  necessarily  imply 
that  it  can  be  obtained  from  the  PY  hard  sphere  density  prof iles . Yet , 
the  general  behaviour  of  the  universal  bridge  functions  for  the  density 
profiles  in  the  range  z>0  is  similar  to  that  of  the  bulk  fluid  and  may 
be  inferred  from  the  solution  of  the  PY  equation  for  hard  spheres  near  a 
wall,  which  is,  quite  naturally,  a  good  starting  point. 

The  solution  of  the  PY  equation  for  a  one  component  hard  sphere 
system  of  bulk  packing  fraction  n,  near  a  hard  wall  yields  [9] 

l  .  H 

-<l+2rtf/(l-rtf  for  z<-l/2 


I  l  i  i  H 

c(z)-  <  [-(l+2r^  +6ry(l+iy'2)  (l+2i^)  (z+1/2)  -2n(l+2n)  (z+3/4)  ] /( 1  -r^ 


for  -1/2<z<1/2 


0 


for  z>l/2 


(6.1) 


where  z  is  is  the  distance  to  the  wall  and  the  hard  core  diameter  d-2R 
is  the  unit  of  length. g(z)  can  be  calculated  from  the  planar  OZ  as 
described  in  [9],  By  a  method  similar  to  that  of  the  bulk  case  [ 1 ] , the 
bridge  functions  for  the  density  profiles  are  given  by 


-c(z)-l-ln[ -c(z)  ;  • 


,z<l/2 


B 


g(z) -l-ln[g(z) ] 


,  z>l/2 


(6.2) 


These  functions  are  continuous , since  as  it  is  well  known 

c(l/2)-g(l/2)-(l+2r\)/(l-il)1  (6.3) 

in  agreement  with  (3.14)  for  hard  spheres 

1/2  g2(l/2)-l/2£(9p/dp) |r  — l/2c(0)|  (6.4) 

The  density  profiles  exhibit  more  pronounced  structure  than  the 

bulk  pair  correlation  functions , but  the  general  shape  of  B  (z)  for  z>0 , 

Tins 

is  very  similar  to  that  of  the  bulk  B  (r)  for  r>0.Note  however,  that 
negative  values  of  g(z)  occur  for  r£>0.45,  while  negative  values  for  the 
bulk  g(r)  occur  only  for  r^0.63  .Bulk  MHNC  calculations,  as  well  as  the 
variational  perturbation  theory  map  the  fluid  range  of  simple  fluids 
onto  [13]  r\_(p , T)< . 45  ,  with  r^(j> , T)-0 . 45  serving  as  a  Lindemann  type 
freezing  rule.  The  singularity  of  the  functions  (6.2),  associated  with 
negative  or  zero  values  of  g(z) .which  limits  their  range  to  r\<0.45,  may 
be  more  than  accidental .Notwithstanding  this  speculative  interpretation, 
one  should  be  cautious  in  applying  the  bridge  functions  (6.2)  for  a 
fluid  near  its  bulk  freezing  point. The  overall  general  features 
regarding  the  application  of  the  MHNC  scheme  to  density  profiles  are 
similar  to  the  bulk  case . Accurate  simulation  data,  when  compared  to  MHNC 


results  will  provide  the  ultimate  test  for  the  accuracy  of  the  bridge 
functions . 


VII -  Cone lus ions 


In  this  paper , devoted  to  the  analysis  of  the  density  profile 
problem  as  the  'planar'  limit  of  a  bulk  mixture  in  which  one  of  the  bulk 
particles  increases  in  size  to  become  the  planar  wall, we  provided  strong 
evidence  for  the  universality  of  the  bridge  functions  for  the  density 
profiles.  Although  we  still  have  room  for  improvement  of  the  parametric 
representation  of  the  bridge  functions,  the  resulting  MHNC  theory  is 
likely  to  be  as  accurate  as  in  the  bulk  case. A  program  for  computing 
density  profiles  by  means  of  the  MHNC  scheme  is  currently  underway. 

An  attempt  to  implement  the  universality  of  the  bridge  functions 
for  bulk  uniform  systems  to  the  treatment  of  nonuniform  fluids  has 
already  been  made  with  considerable  success  [32]. From  the  present 
analysis,  it  seems  however  that  fruitful  information,  leading  to  a 
systematic  improvement  of  the  theory  of  nonumiform  fluids,  may  be 
obtained  by  comparing  the  MHNC  results  for  the  density  profiles  to  those 
obtained  using  the  weighted  density  functional  formalism  [33], 

Final  test  of  the  accuracy  of  this  theory  will  be  achieved  by 
comparison  to  computer  simulations .We  hope  to  perform  such  a  study  in 


the  future . 
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Table  1  :Densi 


profile  sum  rule  for  the  OCP  (see  text,  and  also 
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Figure  caption 

Figure  1  :  The  bulk  bridge  function  parameter  n  as  a  function  of 

F-l/2  [  l+il(dP/dp)  ]  -ftP/jD  for  different  potentials  .  The  data  is  taken  from 

references  [1,13,14],  The  Lennard  -Jones  results  include  data  for  both 


